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Abstract: We compute (as functions of the shape and Wilson-line moduli) the one-loop 
Casimir energy induced by higher-dimensional supergravities compactified from 6D to 4D 
on 2-tori, and on some of their orbifolds. Detailed calculations are given for a 6D scalar 
Qh! field having an arbitrary 6D mass m, and we show how to extend these results to higher-spin 

fields for supersymmetric 6D theories. Particular attention is paid to regularization issues 
| and to the identification of the divergences of the potential, as well as the dependence 

of the result on m, including limits for which m 2 A <C 1 and m 2 A 3> 1 where A is the 
volume of the internal 2 dimensions. Our calculation extends those in the literature to 



very general boundary conditions for fields about the various cycles of these geometries. 
The results have potential applications towards Supersymmetric Large Extra Dimensions 
(SLED) as a theory of the Dark Energy. First, they provide an explicit calculation within 
which to follow the dependence of the result on the mass of the bulk states which travel 
within the loop, and for heavy masses these results bear out the more general analysis of 
the UV-sensitivity obtained using heat-kernel methods. Second, because the potentials we 
find describe the dynamics of the classical flat directions of these compactifications, within 
SLED they would describe the present-day dynamics of the Dark Energy. 
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1. Introduction. 

The Casimir energy for various field theories compactified on 2 internal dimensions has 
been extensively studied (see for example refs. |i[|- [p!q1 ). In this paper we return to this cal- 
culation for the particularly simple examples of 2-tori and their orbifolds. At the technical 
level, our aim in is to provide a generalization of previous calculations for fields compacti- 
fied on T2 and T2/ZN to include very general boundary conditions, including those which 
would be generated by (constant) Wilson-line backgrounds. 
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Our physical motivation for performing this calculation is due to its relevance for 
the recent proposal for using 6D supergravity to shed light on the cosmological-constant 
problem |17], [18| within the context of Supersymmetric Large Extra Dimensions (SLED). 
(See also [19] for related discussions and similar proposals.) Within this proposal the 
extra dimensions must presently be sub-millimeter in size, and the recently-discovered [20] 
cosmological Dark Energy density corresponds to the Casimir energy of the model's bulk 
fields as functions of the moduli of these large 2 internal dimensions. The smallness of the 
Dark Energy density in this picture is ultimately traced to the large size of these extra 
dimensions. In this picture these moduli can be extremely light in a technically natural 
way, and so the Dark Energy phenomenology is controlled by the dynamics of the moduli 
as they respond to this Casimir energy. Although simple cosmologies appear to be possible 
along roughly these lines [21], a more detailed determination of their viability requires a 
more careful calculation of the modulus potential which is provided by the Casimir energy, 
such as we present here. 

Torus-based Casimir energies are particularly well-studied, and recent one-loop studies 
include H, who study the energy of massless scalars compactified on T2 and T2/Z2 at the 
special modular point corresponding to an orthogonal underlying torus. Ref. ||, on the 
other hand, computes the full modulus-dependence of the Casimir energy for various mass- 
less fields compactified on T2, assuming these fields to satisfy periodic boundary conditions 
about the cycles of the background geometry. Ref. computes for T2 compactifications 
the dependence of the Casimir energy on a particular modulus, argt/, but restrict some 
others (by choosing equal toroidal radii). Refs. ||] and |J consider the T2/Z2 orbifold, 
with moduli fixed to those values appropriate for an orthogonal underlying torus. Ref. [10] 
considers the case of T2/Z2 or T2, including the presence of Wilson lines, but only com- 
putes the Wilson-line dependent part of the result. Other recent calculations make similar 
assumptions [11, |l2], 12]. 

In this paper we present results for the fields which appear in supergravity models 
compactified on a 2-torus, T2 and some of its orbifolds, T-z/Zn, as functions of the relevant 
moduli and the higher-dimensional particle mass, to. We do so for a very general set 
of boundary conditions about the cycles of the background geometry, such as could be 
generated by the presence of nontrivial Wilson lines wrapping these cycles. We provide 
formulae which lend themselves to numerical evaluation, and focus in particular on the 
form of the result in the large- and small-m limits. 

A spin-off of this calculation is the information it provides about the large-TO limit, and 
so of the ultraviolet sensitivity of these Casimir-energy calculations. This UV-sensitivity is 
a crucial part of the SLED proposal, and explicit calculations such as those presented here 
provide important checks on the more general, heat-kernel, UV-sensitivity calculations for 
6D backgrounds, presented in a companion paper, These general heat-kernel results 
properly describe the explicit dependence of the toroidal example considered here, and also 
show how flat geometries like tori are dangerous gedanken laboratories, because they are 
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particularly insensitive to large-m UV effects. 

The explicit orbifold calculation also provides the simplest example of the appearance 
of new, brane-localized, UV divergences — a phenomenon which is generic to Casimir- 
energy calculations in the presence of branes due to the singularities and boundaries which 
these branes typically induce in the background geometry [23, 22]. We exhibit this new 
divergent contribution explicitly and show that it arises due to the orbifold projections 
which must be performed. 

So far as they go, our results support the SLED framework inasmuch as all of the 
one- loop contributions we find are at most of order 1/A 2 , where A is the compactification 
area. This makes them no larger than is required for the success of the SLED proposal. 

The paper is organized as follows: In the next section we describe the relevant parts 
of the toroidal geometry and the boundary conditions for fields on this geometry for which 
we compute the Casimir energy. The explicit calculations are first performed in detail 
for complex scalar fields, for which we also display the ultraviolet-divergent and large-m 
dependent parts. Results for the finite part of the Casimir energy are given for general 
boundary conditions for massless scalar fields and these are then extended using a simple 
argument to obtain the corresponding results for massless higher-spin fields in 6 dimensions. 
Section |3| addresses the same issues for the case of T^/Zn orbifolds, with the T2/Z2, T2/Z4. 
orbifolds treated in detail. The small-m/large-m cases are again discussed. The main text 
only quotes the results of the calculations, and full technical details are provided in the 
Appendix. Many of the tools of this Appendix including the detailed calculations of series 
of Kaluza-Klein integrals, can prove useful for other applications such as loop corrections 
to gauge couplings in gauge theories on orbifolds with discrete Wilson lines. 



2. Casimir Energy for 2-Tori. 

In this section we compute the Casimir energy for various 6D fields in a compactification to 
4 dimensions on a 2-torus 7^, for a very broad class of boundary conditions for these fields 
about the two cycles of T2. Although the result is interesting in its own right, it is also the 
starting point for the later calculation of Casimir energy on orbifolds (which exhibit the 
effects for the Casimir energy of the presence of co-dimension 2 branes). 
We define T2 by identifying points on the plane according to 

(2/1,2/2) = (2/1 + n 2 L 2 cos8 + n 1 L 1 ;y 2 + n 2 L 2 sin 6>) , (2.1) 

where ni j2 are integers and 6, L\ and L 2 are the three real moduli of the torus (see Fig. [l]). 
Equivalently, in terms of the complex coordinate z = y\ + iyi this is 

z ^ z + {n 2 U + ni )Lx , (2.2) 

with the complex quantity U defined by U = exp(i#)L 2 /Li = U\ + 1U2, (U2 > 0). 
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Figure 1: 2D torus characterized by the moduli 6, L\ and L 2 , defined by the two identifications 
(2/1,2/2) = (2/1 + £1,2/2) and (2/1,2/2) = (2/1 + L 2 cos9,y 2 + L 2 sm6). 

2.1 Scalar Field Casimir Energy 

Consider a complex 6D field $ on a space-time compactified to 4 dimensions on the above 
2-torus. Writing the six coordinates as {x^,yi}, with p, = 0, ...,3 and i = 1,2, assume the 
scalar satisfies the following boundary conditions 

$(x, vi + n 2 L 2 cos 9 + ni Lf, y 2 + n 2 L 2 sin 0) = e 2ni( - ni pl+n2 ^ §(x, yi, y 2 ) , (2.3) 

with pi,2 being two real quantities. The choices pi j2 = 0, 5 correspond to periodic or anti- 
periodic boundary conditions along the torus' two cycles. More general values of pi are also 
possible, such as when $ transforms non-trivially under a gauge group for which nonzero 
Wilson lines are turned on for the corresponding cycles. 1 Here we consider pi >2 arbitrary. 

Expanding the scalar field in terms of eigenfunctions of the 2D Laplacian, D 2 = d 2 +0 2 , 
according to 

y 2 ;pi,p 2 ), (2.4) 

ni,n 2 

we have □ 2 /m,n 2 = - M n l5 n 2 /m.rw with 

^n 1 ,n 2 (Pi,P2) = ^-\n 2 + p 2 -U(n 1 + Pl )\ 2 . (2.5) 
Here *4 = LiL 2 sin# denotes the area of the torus. The mode functions are given by 

f ni ,n 2 {yi,y 2 ;pl,p 2 ) = e 2 ^[(™i+^)(^-^cote)/L 1 +(n 2 +p 2 )j /2 /(L 2 sme)] _ ^ 



To see this (for an orthogonal torus) notice that the toroidal boundary conditions preclude removing 
a constant gauge potential, such as A\ = a, using only strictly periodic gauge transformations. (This 
corresponds to the Wilson line W = f^ 1 A m dy m — aL\.) However, A\ can be removed using a singular 
gauge transformation having parameter u) = ayi, at the expense of changing the boundary conditions of 
charged fields: $(2/1 +L%) = e lqaLl <3>(j/i ) , where q is $'s charge. We see from this that 2npi = qaL\ — qW. 
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We now compute the vacuum energy for a complex scalar field having 6D mass m, com- 
pactified on T2 with the above boundary conditions. Denoting the vacuum-energy per unit 
3- volume for such a scalar field by V(pi,p2), we have 



d d p 
-, — rr In 

27T 



d 



f> 2 



ni )2 eZ ' 

where we have continued the momentum integration to Euclidean signature, and we reg- 
ulate the ultraviolet divergences which arise in the sum and integral using dimensional 
regularization, with the complex quantity 2 d = 4 — e ultimately being taken to 4. Here 
\i denotes the arbitrary mass scale which arises in dimensional regularization, and which 
drops out of all physical quantities. 

A potential subtlety arises in the above expression for massless 6D fields (m = 0) if 
both p\ or P 2 are integers, because in this case there is a choice of integers (ni, n-z) for which 
M niTl2 vanishes. In this case it is convenient to keep m nonzero so that all manipulations 
remain well-defined, with m taken to zero at the end of the calculation. In principle, 
one must be alive to the possibility of unexpected singularities appearing when m tends 
to zero after renormalization, such as the familiar infrared mass singularities of Quantum 
Electrodynamics ]24[| . As usual, such infrared problems are less severe in higher dimensions 
and we shall see that there is no such obstruction to taking m — > for our applications in 
6 dimensions. 

The calculation of the two infinite sums in ( |2.7D is tedious, and is given in detail in 
Appendices [A] and g — c.f. eqs.([A^lD, (fETID , ( |B-14[ ) to QB-18D — using the approach of 



|25| ] . In what follows we quote only the final results which are appropriate to the discussion 
at hand. The next three sections respectively concentrate on the ultraviolet-divergent part 
of the result, as well as the finite part in the cases where the 6D scalars are either massless 
(m 2 A — > 0) or very massive (m 2 A 3> 1). 

2.1.1 Ultraviolet Divergences for 2-Tori 

The ultraviolet divergent part of V in (|2.7|) denoted Vqo is, with e = 4 — d (see eqs. ( |B-l| ), 



( grip to ( pnsD ) 



V 00 {p liP2 ) = — rr (2-f 



which is valid for arbitrary m. (Eq. (|2.8| ) shows the importance of keeping a non-zero m 
when discussing ultraviolet divergences in dimensional regularization (DR), since these can 
easily be missed if m = 0.) This expression has several features on which we now remark 
(and which agree with the more general analysis of the ultraviolet divergences in 6D field 
theories compactified on Ricci-flat backgrounds given in a companion paper p2|). 



2 The conventions here differ from those used in our previous ref. |22j where d = 4 — 2 e. 
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• First, the divergent part depends on the moduli, an d 0, only through the toroidal 
area A = L\Li sin#, and is interpreted as being a renormalization of the 6D cosmological 
constant. Note also that the UV divergence vanishes if we take 9 — > 0, corresponding to 
collapsing the two cycles of the torus onto each other down to one dimension less. 3 This 
agrees with the well-known absence of one-loop UV divergences for a broad class of theories 
when they are dimensionally regularized in odd dimensions. 

• The proportionality to m 6 is also what is required on dimensional grounds (in dimensional 
regularization) for a contribution to the 6D cosmological constant. The absence of other 
powers of m, such as an ^.-independent result proportional to m 4 , is a consequence of the 
torus being flat, and is not true for more general curved spacetimes [p2| . Once the UV 
divergence is renormalized into the 6D cosmological constant there is no obstruction to 
taking m — > 0, unlike the situation for massless 4D theories. We consequently feel free to 
simply set m = in subsequent applications of our formulae for V mn where 

yren _ y _ y^ (39) 

with V as in ( |2.7D and U^ its divergent part. 

• The divergent part of V is pj-independent and so does not depend on the boundary 
conditions of the scalar field, eqs.([2.3|). Again this agrees with general arguments, since the 
short-wavelength modes responsible for the UV properties are not sensitive to the boundary 
conditions which depend on the global properties of the background geometry. We shall 
see that for orbifolds new divergences are present corresponding to counterterms localized 
at the fixed points, and these new divergences can depend on the nature of the boundary 
conditions imposed on the covering space. 

• The 1/e pole which appears here represents a bona fide 6D divergence. This is at first 
sight surprising, since e represents the difference between d and 4 rather than 6, and it is 
introduced for each of the Kaluza-Klein (KK) modes. To understand this it is important 
to recognize that our expressions contain two separate sources of UV divergence: the 
integration over 4-momentum, p, and the two sums over KK mode numbers, n^. In our 
calculations the dimensional continuation is e = 4 — d away from four in order to regularize 
the p-integration. On the other hand, the KK mode sums are managed using zeta-function 
techniques, and the presence of e ensures the regularisation of these sums. With these 
choices the leading inverse powers of e obtained turn out to be precisely those which would 
be obtained starting from 6D and following the powers of e' = (6 — d) . This equivalence is 
shown in more detail for an explicit example (using a spherical geometry, for which more 
divergences may be followed) in ref. |p2| . 

We now examine the finite parts of the Casimir energy density. 

3 This limit must be treated with care, however, since the calculations of Appendix ^ also require 
U2 ~ sin# to be finite and nonzero. 
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2.1.2 Massless Fields in 6D 



In the massless limit, m — > 0, the result for the vacuum energy density is (see eqs,( A-3| ), 

@, (EH to (El)) 



+ E 

nieZ 



m=0 



\4 2 I 90 



2i - A P, ( 1 -5A^ 1 -2A^ i +6A3 i ) 



(n 1 + A pi ) 2 Li 3 K 1 ) + 3|ni + A ^ 1 



2vr^ 2 



Li 4 (o- ni ) + 



(2vr?7 2 ) : 



Li 5 (cr ni ) + c.c. 



with 



Cni>0 



2i tt( A P2 -UA P1 -C/ni ) 



and 



1 



Cni<0 — — 



(2.10) 



ff ni>0 

where a ni denotes the complex conjugate of a ni . In these expressions < A Pi < 1 
represents the fractional part of p«, as in pi = [pi] + A Pi , where [pi] G Z is the largest 
integer smaller than or equal to pi. The poly-logarithm functions which appear here are 
defined by the sums p6| 

X r ' 

n° 



Li CT (x) = ^;^. 



(2.11) 



n=l 



The point of rewriting the initial two sums into the ones written here is that these converge 
well and so are useful for numerical purposes. Figure |2| plots V(pi,P2) as functions of the 
moduli Ui and U2 for various choices for the boundary conditions (pi,^)- 

We have checked that the above formula agrees with the particular cases studied in 
the literature. For instance in the special case p\ = P2 = we find 

v z rii=l 



+^ Li 4 (^) + Li 5 (g ni ) + cc. 



2irU 2 



4tt 2 U 2 



(2.12) 



where q = e 2l7TU . This agrees with the result given in ref. 



2.1.3 Heavy-Mass Dependence 

The generality of the calculation in Appendix |A|, |b| also allows the explicit exhibition of 
the heavy-mass limit, m 2 A —> 00, of the Casimir energy. This is of particular interest for 
Supersymmetric Large Extra Dimensions, where the naturalness of the description of the 
Dark Energy density relies on the Casimir energy only depending weakly on the masses of 
heavy fields in the 6D bulk. Using formulae ( p3-15|) , ( B-16| ), ( B-18| ) of the appendix it may 
be shown that if m » {1/Li, 1/(1/2 sin 6*)}, leading to m 2 A » 1, then 



V(pl,p2) 



m°A /2 
384vr 3 \1 



In 



(2.13) 
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thus powers of m 2 A other than m 6 are exponentially suppressed. These expressions agree 



well with the general results of ref. [22], which identify the large-m behaviour using general 
heat-kernel techniques. For general geometries there can be powers of m in the large-m 
limit, but the leading such powers are proportional to local effective interactions which 
involve polynomials of the background fields and their derivatives. For the simple toroidal 
geometries considered here all of these local interactions vanish, leading to the exponential 
mass suppression found above. 

The absence of powers like m 4 or m 2 in the large-m limit is more difficult to understand 
from the point of view where the 6D calculation is regarded as simply being the sum over an 
infinite number of 4D contributions, each of which can themselves have such powers of m. 
As is clear from the general 6D analysis of |23] , the absence of these terms may be traced to 
the requirements of locality and general covariance in 6 dimensions — requirements which 
are easily missed in a KK mode sum calculation. 

2.2 Higher-Spin Fields on T2. 

With an eye towards applications to supersymmetric theories, in this section we compute 
the corresponding results for the Casimir energy for other massless fields in 6 dimensions. 
We do so using the trick of ref. H , which uses the prior knowledge that the Casimir energy 
must vanish once summed over the field content of a 6D supermultiplet, provided that 
these fields all share the same boundary conditions about the cycles of T2 (and so do not 
break any of the supersymmetries). 

To this end we reproduce as Table |l] a table 
from ref. |22| listing the massless field content 
of some of the representations of (2, 0) super- 
symmetry in 6 dimensions. In this table the 
scalars are real, the spinors are symplectic-Weyl 
and the 2-form gauge potentials are self-dual or 
anti-self-dual. 

The argument of ref. Q uses the observa- 
tion that a single symplectic-Weyl fermion and 
two real scalars preserve 6-dimensional (2, 0) 
supersymmetry in a toroidal compactification 
for which they share the same boundary condi- 
tions, (pi,p2)) about the torus' two cycles. The 
Casimir energy for these fields must therefore 
cancel in order to give a vanishing result for 
the contribution of a hypermultiplet. Since we 
know the scalar result for general pi, we may 
infer from this that the Casimir energy for a 
single symplectic-Weyl fermion must be precisely —1 times the result quoted above for a 
complex scalar field having the same boundary conditions. 



Multiplet 


Field Content 


Hyper 


W>+, 20) 


Gauge 


(A M , 2^_) 


Tensor 


(A+ N , 20+, 0) 


Gravitino I 


[i>M+, 2A Mn) 


Gravitino II 


(ipM~, 2A M ,4>-) 


Graviton 


(9m N, 2-0 m+, A MN ) 



Table 1: The field content of some of the 
massless representations of (2,0) supersym- 
metry in 6 dimensions. The fermions are 
taken to be symplectic-Weyl, and their plus 
and minus subscripts correspond to their chi- 
rality, as measured by their I\ eigenvalue. 
The + (— ) sign on the 2-form potential simi- 
larly indicates it is self (anti-self) dual. There 
is also an equivalent set of representations 
where all plus and minus signs are inter- 
changed. 
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Using the identical argument based on the vanishing of the Casimir energy summed 
over the field content of a gauge multiplet similarly shows that the Casimir energy of a 
6D gauge boson must be —2 times the result for a 6D symplectic-Weyl fermion — and so 
is +2 times the result for a 6D complex scalar — having the same boundary conditions. 
Arguing in this way allows the inference of the Casimir energy for all of the other fields 
appearing in the supermultiplets listed in Table |]. The results found in this way are 

Vi/2{P\,P2) = -V{pi,p 2 ) 

Vi( Pl ,p 2 ) = 2V( Pl ,p 2 ) 
3 

V K r(pi,p 2 ) = ^V(pi,p 2 ) 
V3/ 2 (px,p 2 ) = -3V(p 1 ,p 2 ) 
V 2 (px,p 2 ) = -V( Pl ,p 2 ), (2.14) 

where the divergent and finite parts of the right-hand side of this equation are given 
explicitly by eqs. ( |2.8| ) and (2.10), above. Here Vi/ 2 , V\, Vkr, V%/ 2 and V 2 respectively 



denote the results for a symplectic-Weyl fermion, a gauge boson, a Kalb-Ramond self-dual 
(or anti-self-dual) 2-form gauge potential, a symplectic-Weyl gravitino and a graviton. 
Given these expressions, it is simple to compute the nonzero Casimir energy which re- 



sults when 6D supersymmetry is broken d la Scherk and Schwarz [27], by assigning different 
boundary conditions to different fields within a single supermultiplet. Supersymmetry is 
broken in this case by the boundary conditions themselves. For instance, if the symplectic- 
Weyl fermion in a hypermultiplet has boundary condition (pi, p 2 ) but the complex scalar 
has boundary condition (pi, p 2 ) then the Casimir energy for this hypermultiplet would be 

V hyper = V 1/2 ( Pl ,p 2 ) + V(p[,p' 2 ) = V{p' x ,p' 2 ) - V{ PuP2 ) . (2.15) 

Similarly applying different boundary conditions to the constituents of a gauge or tensor 
multiplet gives 

Ugauge = V 1 ( P1 , P2 ) + 2 V 1/2 (p[, P ' 2 ) = 2 [V( Pl ,p 2 ) - V(p[, P ' 2 )\ 
U tens0 r = V KR ( Pl ,p 2 ) + 2 V 1/2 (p[ , p' 2 ) + - V(p'(, p' 2 ') 

= \v{p 1 ,p 2 )-2V{p' 1 ,p' 2 ) + \v{p>lp'l) , (2.16) 

and so on. 

Figure || gives in addition to the plots of V(p\, p 2 ), the differences V(p\, p 2 ) — V(p[, p' 2 ) 
as functions of U\ and U 2 for various choices for the boundary conditions (pi,p 2 ) and 
(p'i,p 2 ). The periodicity wrt U\ in the plots is a remnant of the SL(2,Z)u symmetry 
(modified by non-zero Wilson lines). For U 2 = L 2 /L\ sin# > 0(1), one has flat directions 
for V (as function of U\,U 2 ). This changes for U 2 <^ 1 (say if <C 1) when V develops 
maxima/minima. For values of pi other than those in the figure, the peaks in these plots 
have different height. 
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(a) The potential V ren for p x = p 2 = 0. 



(b) The potential V ren for p x = 1/2, p 2 = 0. 




Figure 2: The (finite part of the) potential V, up to an overall positive factor and for fixed area 
A. The plots are in function of U\, U2, (U = U\ + HJ2 = L^/Lx exp(i#)) for various p\ : 2- See also 
eqs.( ^.10| ) and ( 3.18 ). The peaks of the plots indicate (moduli) divergences present at U2 <C 1 when 
the two dimensions collapse onto each other (8 = 0). 
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3. Casimir Energy for Orbifolds 



None of the previous results included the effects of 3-branes within the bulk when computing 
the Casimir energies. We now extend these calculations to some simple examples which 
include branes, and for which the background geometry includes the back-reaction of the 
brane tensions by incorporating the appropriate conical singularities at the brane positions. 
We only consider here brane singularities which correspond to the specific defect angles 
which arise when an orbifold is constructed from the 2-torus by identifying points under 
the action, Zn, of a discrete set of rotations. We analyze separately the case of T2/Z2, with 
full details, and then the orbifold T2/Z4 whose technical details differ considerably. For 
T2/Z4 we use a general method which can be applied to the remaining T2/Z3 and T2/ZQ. 

To this end we return to the description of 72 as a complex plane, z = y\ + iy2, 
identified under the action of a lattice of discrete translations as in eq. (pi!). Following 



standard practice, we construct an orbifold from this torus by further identifying points 
under the action of the Zn rotations defined by 

z^r k z, (3.1) 

where r = e 27ri / N . This gives a well-defined coset space, O = T^/Z^, provided that these 
rotations take the initial lattice which defines the torus onto itself. Notice that if JV = 2, 
then the rotation z = — z is automatically a symmetry of the lattice for any value of the 
moduli L\, L2 and 9, and so these three quantities are also moduli of the resulting orbifold, 
O. On the other hand, \i N > 2 then the rotation is a symmetry of the lattice only for 
specific choices for the complex structure: U = {L2/Li)e td = r, and so L2 = L\ = L and 
6 = 2ir/N, and so only one modulus, L, in this case survives. 

The coset T^/Z^ is an orbifold rather than a manifold because of the metric singu- 
larities which arise at the fixed points of the group. For instance, in the case T2/Z2 there 
are 4 such points, corresponding to z = 0, ^, ^ U and ^(1 + U). The metric has a conical 
singularity at each of these points, whose defect angle is n. For further details on orbifolds 
and their fixed points see Appendix IDL 

3.1 A Scalar Field on T 2 /Z 2 . 

We now compute the Casimir energy for a complex scalar field, compactified from 6D 
to 4D on the orbifold O = T2/Z2. As before, the 6 coordinates are taken to be {x, y{\, with 
the orbifold corresponding to the coordinates yi,i = 1,2. We consider the 6D scalar field 
to satisfy the boundary conditions 

$(x, yx + Li; y 2 ) = e 2mpl $(x, y u y 2 ); pi = 0, 1/2. (3.2) 
$(x, Vl + L 2 cos 6; y 2 + L 2 sin 6) = e 2mp2 <f>(x, yi ,y 2 ); p 2 = 0, 1/2. (3.3) 
$(s, -yi, -y 2 ) = ± $(x, yi,y 2 ) . (3.4) 



Condition (|3.4| ) is possible because of the new cycle that the orbifold has (which the torus 
does not). As is indicated in eqs. ( |3^ ) and (|Q|), the quantities pi are no longer free to take 
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any real value in this case, because the underlying Wilson lines must be compatible with 
the orbifold rotation. For instance, when acting on the coordinates it is straightforward 
to show that the composite transformation I = SoPoS gives X(?/i, 2/2) = (— 2/1, —2/2), if 
^ (2/1,2/2) = (2/1 + £1,2/2) and P(2/i, 2/2) = (—2/1,-2/2)- Applying the same transformations to 
and using the boundary conditions fl3,2|) and G3.4| ), consistency requires exp(4-7ripi) = 1, 
and so pi = re/2 for some integer re. A similar argument implies that P2 is also half-integer. 
If we require, without loss of generality, < pi < 1, then we see that consistency requires 
p\ 2 = 0, 1/2. Altogether there are 8 possible choices for the boundary conditions, denoted 

by (pi,P2) ± - 

Because the orbifold reflection is a symmetry of O2, these reflections have a natural 
action on its eigenfunctions, f ni ,n 2 - Using the explicit expressions obtained earlier, eq. (^), 
for the toroidal mode functions 

VA 

= fn 1 ,n 2 (yi>y2;pi,P2) , 

= fn' 1 ,n! 2 (yi,y2]Pl,P2) , (3.5) 

where re- = — rii — 2pi ensures ■n! i + pi = — (rij + pi). Notice that ■n! i defined in this way 
remains an integer because for the Z2 orbifold pi = 0, |. 

Using this action it is straightforward to specialize the toroidal mode expansion, 
eq. (|2.4j ), to fields on T2/Z2 with boundary conditions (p\,p2) ± - We now write these 
expansions explicitly, in terms of the real and imaginary parts of the mode functions 

(1). (pi,p 2 ) = (0,0) ± . In this case f nim (-y\, -1/2) = ,-„ 2 (yi, 2/2) and so using the 
boundary conditions of eqs. (|3.4[) in the mode expansion of eq. (|2.4|) gives 



2 

^ >+ (x,y 1 ,y 2 ) = V -g— -4>„ u n 2 (x) e nin2 (y 1 ,y 2 }0,0), 

ni>0, n 2 GZ 

2i 

* (a;, 2/1, 2/2) = tx— :<j>mm( x ) 9m, n 2 (2/i,2/2; 0,0), (3.6) 

m>0, 7i2£Z 

where the superscript on $ indicates the sign chosen for the orbifold projection, and 5 nit o 
is the usual Kronecker delta-function which vanishes unless m = 0, in which case it equals 
unity. Notice that because 2/0,0(2/1, 2/2; 0, 0) = 0, the mode (ni,re2) = (0,0) is absent in the 
double sum for <E>~, as is indicated by the primed double sum. 

(2). (pi,p 2 ) = (0, l/2) ± . In this case / ni ,n 2 (-2/i, ~V2) = f-m,-n 3 -l(yi, 2/2) and so the 
mode expansion becomes 

2 

^ + (x,yi,y 2 ) = V 7^— -^ni.na^) e ni ,„ 2 (2/i,2/2;0, 1/2), 

ni>0, n 2 eZ 

2i 

$~(x, 2/1, 2/2) = V -t— -0 m , n2 (ar) (7 ni n2 (2/1, 2/2; 0,1/2). (3.7) 

ni>0,n 2 eZ 
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(3). (pi,p 2 ) = (1/2, O)*. Here f ni ,n 2 (-yi, ~V2) = /-m-i -n 2 (2/i, 2/2), and so 



(x) {yuM 1/2,0), 

ni>0, «2eZ 

^ _ (x,2/i,y 2 ) = ^ 2i(f) nun2 (x) g ni ,n 2 (yi,y2; 1/2,0), 

ni>0, n2GZ 



(3.8) 



(4). (pi,p 2 ) = (1/2, l/2) ± . In this case f ni ,n 2 (-yi, -2/2) = /-m-i ,-n 2 -i(yi, 2/2), and so 



Vni,n2 

(X) °ni,n2 

(j/!,^; 1/2, 1/2), 

ni>0,n 2 eZ 

^ ( x ,yi,y2) = ^2 2i0 niin2 (x) <7ni,n 2 (2/1,2/2; 1/2, 1/2). 
ni>0, «2€Z 



(3.9) 



There are two ways to compute the Casimir energy for the scalar field $ on T2/Z2 with 
these boundary conditions. One approach is to recognize that the scalar propagator on the 
orbifold may be obtained from the propagator on the torus using the method of images: 

G%{x -x',y- y') = G r {x -x',y- y') ± G r (x -x',y + y') , (3.10) 

and following the implications of this for the vacuum energy. The second approach is to 
directly perform the KK mode sum over the modified mode functions given above. Both 
lead to the same result, and we present the mode-function derivation here because, albeit 
more involved, it can be extended to the case of T 2 /Zn and allows a general discussion of 
the ultraviolet divergences for T 2 /Zn. 

The vacuum energy density per unit 3-volume written as a mode sum is given by 



Vb(0, 0) ± 
Vb(0,l/2) ± 
Vfe(l/2,0) ± 
Fo(l/2,l/2) ± 



(-)' ^4-d r 

Jr* C 7 2*"i-° J 

ni>0,n2G^ 



E 



ni>0, ri2GZ 



d d p 
(2vr) d 

d d p 

2 5 "i.° J (27r) d 



4-d 



In 



In 



p 2 + MZ un2 (0,0)+m 2 
M 2 



P 2 + M^ lin2 (0,l/2)+m 2 



m>0,n 2 eZ 

E " 

ni>0, ri2SZ 



4-d 



d d p [p 2 + M2 iin2 (l/2,0)+m 2 



4-d 



(2vr) rf 
d d p 

{2ll) d 



In 



In 



p 2 + M ra 2 1;n2 (l/2,l/2)+m 2 
/, 2 



, (3.11) 



where d = (4 — e) and the symbol (— )' on the sum in the first line indicates the exclusion 
from the sum of the single mode (ni, n 2 ) = (0, 0), but only for the case of (0, 0)~ boundary 
conditions. 

As might be expected from the approach based on the method of images, these expres- 
sions may be evaluated in terms of the corresponding quantities on the torus. To see this for 
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the mode sums we denote the summands of these expressions by W(|ri2 + P2 — U (n\ + pi)| 2 ) 
in order to emphasize their dependence on mode numbers and moduli. It is then simple to 
use the invariance of W under changes in sign of (rtj + pi) to prove the following identities: 



Hi 2 



E ^ W (K"^i| 2 ) =\ E ^(K-t/n 

;0;ri2£Z 7ii,n2eZ 

E + =\ + 1/2 -Un 



2 

ni>0;n2SZ ni,ri2SZ 



2 

II 



w(\n 2 + p 2 -U{n 1 + l/2)\ 2 ) = ~ £ w(|n 2 + p 2 - ^7(m + 1/2)| 2 ) (3.12) 

ni>0;n,2GZ m,ri26Z 

where /)2 is either or ^ in the last line. These expressions allow the derivation of the 
following expressions for the Casimir energies in terms of the toroidal results, V(p\,p 2 ) 

vb(o,o)- = \ [v(o,o)-v zl 

Vo(pi,P2) ± = 2 V(pi,p2) for all others. (3.13) 

Here V zm is the contribution to the torus Casimir energy of the "zero mode" (m, n 2 ) = (0, 0) 
with (pi,p 2 )~ = (0,0)~ (for its expression see ( B-17p and ( B-18|) ). We can now present in 
detail the divergent and finite parts of the sums and integrals in ( |3.11| ), Q3.13 ) . 

3.1.1 Ultraviolet Divergences 

We isolate the divergent part of Vo in eqs.( |3,li| ) and write 

Vo(pi,P2) ± = V ,oo(pi,P2t + V e »(p 1 ,p 2 ) ± (3.14) 

where all divergent terms are included in Vo,oo- As is clear from eqs. ( |3.13| ), for all choices of 
boundary condition except (0,0) ~ on the orbifold the ultraviolet divergences encountered 



are precisely half of those encountered on the torus, eq. (2J 



^^ 2)±= 3^ = T^' if (^2)^(0,0)-. (3.15) 

This divergence may be absorbed, as usual, into a renormalization of the bulk cosmological 
constant. Notice that its coefficient is the same as was obtained earlier for the torus, once 
the divergence is expressed in terms of the area of the orbifold, Ac, which is half the area, 
.A, of the covering torus. 

By contrast, the exclusion of the zero mode for the specific choice (0,0)~ introduces a 
new type of divergence which was not encountered for the torus. In this case the orbifold 
and toroidal divergences differ by the contribution of the n\ = n 2 = mode alone, and 
thus, using eqs.flA^), Q), dB^l4| ), ( p^T5|) , one has 
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The presence of the last term is consistent with the general heat-kernel analysis. Because it 
is proportional to m 4 and is independent of the bulk moduli, it has the right properties to 
be interpreted as a renormalization of the tension of the branes whose presence at the fixed 
points is responsible for the conical singularities in the bulk geometry at these points. As 
before, the ultraviolet-finite part of the Casimir energy obtained after this renormalization 
is nonsingular in the m — > limit, and so we are free to take this limit explicitly in the 
renormalized result for massless 6D fields which we quote below. 

It should be emphasized that although this divergence renormalizes the local brane 
tensions, it arises due to the functional integration over bulk fields. This is a feature which 
arises quite generically for quantum effects in the presence of boundaries and defects, 
whose origin can be understood in detail as follows. The bulk vacuum energy density, 
tMN = {Tmn) is ultraviolet finite (for the flat orbifold under discussion) after the bulk 
cosmological constant is appropriately renormalized. However although tMN is finite, it 
is also position-dependent due to the presence of the orbifold singularities breaking the 
translation invariance of the underlying torus. In particular, £mjv(z/) typically goes to 
infinity as the singular points are approached in a way which diverges once integrated 
over the volume of the orbifold. It is this new divergence which is renormalized by the 
brane-tension counter-term localized at the singularity. 



less complex 6D scalar field compactified on T2/Z2, for the various boundary conditions 
(pitPz) • It is noteworthy that V™(0, 0) vanishes as m — > in dimensional regularization, 
and so the orbifold result is half of the appropriate toroidal result for all choices of bound- 
ary conditions. After the renormalization of the ultraviolet divergences described above, 
one has 



3.1.2 6D Massless Fields 



Using eqs.( 3.12j ) 



we can now give the explicit results for the Casimir energy of a mass- 





(3.17) 



and that 






r "i 

+ ^[n?Li 3 (-^) + ^Li 4 (-^) + 



(2tt[/ 2 ) 2 



3 



Li 5 (-g ni )+c.c. 



ni=l 



- 15 - 



V X 7 711=0 

+ 3(n !t 1/2) Li 4 (^ +1 / 2 ) + t^-L^ Li 5 (^+V2) + c . c ; 



2vr?7 2 



(2vr[/ 2 



l£ en (l/2,l/2) 



m=0 



1 



31 \ (2nU 2 



672 / 180 



oo 

+ £ [(m + i/2) 2 Li3(-g 



m+l/2\ 



ni=0 



3(m + 1/2) _ gm+1/2) + _3 ^ _ gm+ l/2 )+c _ c _ . (3 _ 18) 



2irU 2 



(2ttU 2 Y 



where q = e and the complex conjugate applies only to the series of polylogarithms. 

3.1.3 Heavy-Mass Dependence 

The divergences of the Casimir energy for large m are identical to those for the case of small 
m discussed in Section |3.1.1| . Further, because the orbifold results are simply expressed 
in terms of the toroidal ones, the heavy-mass dependence of the toroidal expressions carry 
over immediately to the orbifold Casimir energy. In particular, in dimensional regulariza- 
tion (and after modified minimal subtraction) the finite parts of the Casimir energy fall 
exponentially for large m, and the only strong m-dependence arises in the divergent terms, 
including the new m 4 term which arises for some of the boundary conditions. 

3.1.4 Higher-Spin Fields on T 2 /Z 2 

The results for massless higher-spin fields on the T 2 /Z 2 orbifold can be read from their 



toroidal counterparts of Section 2.2. This is possible because the orbifold identification 
does not break supersymmetry provided that all of the fields within a 6D supermultiplet 
satisfy the same boundary conditions. The results for the Casimir energy of higher-spin 
fields may therefore simply be read off by multiplying the expressions by the factors 

given in eqs. (2.14). 



3.2 The Orbifold T 2 /Z N with N > 2. 

In this section we outline the steps for computing the Casimir energy for a complex scalar 
field compactified on the T 2 /Z^ orbifolds, with N > 2. Recall that for these orbifolds 

U = e 2i^/iv L 2 = Li = L, 9 = 2it/N (3.19) 

We take the following action of the translation and orbifold Z^ symmetries on the field $ 

<S> 9 (x,T k z) = g k <S> 9 (x,z), 
$ s (x,z + Li) = e 2i7Tpi $ 9 (x,z) 
<£> 9 (x, z + UL ± ) = e 2inp2 $ 9 (x, z) , (3.20) 
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where we use complex coordinates z = y\ + iyi and as before r = e 2ln ^ N and k = 0, N — 1. 
Here 5 is a particular representation of the Z^r transformation, acting on The super- 
script 'g' emphasizes that this representation is not unique, and the explicit form taken by 
<I> S in general depends on which g is chosen. As in the case of T2/Z2, this realization only 
faithfully reproduces the symmetry for specific choices for the pi, whose values we now 
determine. 

The consistency conditions for the action on <I> 3 are found by combining the above 
expressions and using geometrical relations which state how some of the rotations 
can also be expressed as translations on the covering torus. To display these we use the 
complex coordinate z = y% + iy2, in terms of which the lattice of translations which defines 
the underlying torus is generated by ei = L\ and e2 = UL\. Then, depending on the group 
of rotations, Zn, which is of interest, the following restrictions can arise. 

• If an orbifold rotation takes e\ to ei — i-e. there is an integer < k < N for which 
T k e\ = e2 (or r k = U) — then using eqs. ( 3.2CQ to evaluate Q 9 (x,T k (z + ei)) = 
<5 9 (x, r k z + €2) implies 

g k exp(27npi ) = g k exp(2mp 2 ) , (3.21) 
and so we may take p\ = P2 without loss of generality. 

• If an orbifold rotation takes — > — then a similar argument implies 

g k exp(2vri / 9 i ) = g k exp(-2vri / 9 i ) , (3.22) 
and so we may take pi = or 1/2. 

For instance, only the second of these conditions applied to the Z2 orbifold considered 
previously. By contrast, both conditions apply to the case of tt/N rotations which give 
the orbifold T2/Z2N, and so for this case we must take p\ = p% = 0, |. For the Z3 case, 
on the other hand, the first condition applies but instead of the second condition one has 
e2 + T k e2 = — ei, and so we find P2 = pi = 0, 3 or |. These results express the quantization 
on these orbifolds of the underlying Wilson lines which are responsible for the boundary 
conditions which are expressed by the pi. For a more detailed description of Wilson lines 
and their values on orbifolds see Appendix [D]. 

To determine the action of the symmetries ( |3.20 ) on the toroidal mode functions, we 



adapt the discussion of Section 2.3 of ref. [28[] to include the general phases px t 2- It is 
convenient for these purposes to rewrite eq. (|2.6D in complex coordinates 



f n n (~~) = 1 c \(n -r h U)z-(n -r h U)z]/(2L 1 U 2 ) ^ (3.23) 



where = rii + pi, for i = 1,2. The construction of the mode functions for the orbifold 



T2/ZN is done by observing that all of the arguments in [28 remain valid if is replaced 
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by nj = Ui + pi. The basis functions one is led to are given by 

JV-l 



hi 



V fc=0 



(3.24) 



where r/ nj o = 1/y/N if rti = H2 = 0> an d otherwise equals 1. Using this definition, one can 
check that the mode functions satisfy the orbifold condition 

h 9 n 1 ,n 2 (r k z)=g k hl^ 2 {z). (3.25) 

As outlined in [^] not all the functions hn n are independent, since they are related by 

K^ )< ^)=K 1 ,r k y*), (3-26) 

where oj k (r) • (n.i,n 2 ) is defined as a rotation which takes —n x + rn 2 into — + rn 2 = 
21i + r Il2)- The set of all such rotations (i.e. for k = 0,1,.., N — 1) identify N 
domains whose union covers the whole complex plane (Ox; Oy) defined by (Ox; Oy) = 
(— n 1 + T\n 2 ;T 2 n 2 ). Each such domain fixes the set of levels n_i,n 2 which identify the 
set of independent /in 15 n 2 which are not related by the rotation u of ( 3.26| ). This gives 
Rl =ni+pi < and n 2 = n 2 + P2 > as an independent set. Since here < pip < 1, one 
concludes that the conditions n\ < 0, n 2 >0 define an independent set of functions hn^n^ 
for the orbifold T 2 /Z^. 

Using the above considerations, one has the following mode decomposition for 

<f> 9 (x,z) = J2 ®k,n 2 ( x ) h n 1 ,n 2 (z) + -= ^ Hi = ^ + P i , (3.27) 

ni<0,n 2 >0 ^ -A 

which satisfies the desired condition & 9 (x, r k z) = g k z). For example, for T 2 /Z% 

2 



g = T 



g = T 



g = T 



o 



ni<0,n2>0 



<f°(x,z)= Yl ^LnM 4fEw^) 

ni<0,n 2 >0 LVo fc=Q 

r i 2 

fc^x.z) = ^ 7f^ T kf ^( T 



+ 



1 



<o(*) 



rti<0,n 2 >0 LVo fc=o 



fc=0 
2 



(3.28) 



The main difference from the T 2 /Z 2 orbifold is that for T 2 /Z^ the sum over the Kaluza- 
Klein levels is restricted to positive/negative values of ni,2, unlike in T 2 /Z 2 where one sum 
could be extended to the whole set Z of integers. The above mode expansion leads to the 
following expression for the Casimir energy of a complex 6D scalar field on T 2 /Zn, N>2. 



Vo [pi, Pi) = Yl )' 

ni<0, rc,2>0 



4-d 



d d p 
(2n) d 



hi 



P 2 + M n u nM>P2)+ m2 
p? 



(3.29) 



- 18 - 



where for each g one uses the values of pi t 2 which respect the consistency conditions. 

The above domain of summation for ri\ 2 makes the analytical calculation of Vo(pi, P2) 
more difficult than in the case of T2/Z2. The difficulty is caused by the fact that none of 
the sums over n\ , 712 can be extended 4 to a sum over the whole set Z of integers (as we had 
for T2/Z2, eq.( 3.11| )). As a result no (Poisson) resummation of individual contributions to 
V(pi:P2) is possible and the calculation is then more tedious. Although one may still be 
able to work on the covering torus 5 rather than in the orbifold basis, the approach below 
(being valid for any pi^) allows a simultaneous analysis of all orbifolds T2/ZN, N>2. 

After a long calculation (see Appendix |C], eqs.( p-l| ) to ( |C-9| )) one has for V(pi,P2) of 
(3.29) the following result (which is valid for m -C 1/L) 



Vo(pi,p2] 



2 A 2 



V + 



m 6 A 3 -2 
768 vr 3 ~T 



(27TU 2 f 

180 



__ p 2 (1 _ 5p 2_ 2p 4 + 6/9 3 ) 



+ E 

ni<0 



/ \2t- / \ ^ \rii + pi . . . 3 T . / \ 
(ni + pi) L13 (cr ni )H r— p7 L14 {a m ) + 2 L15 (a m ) + c.c. 



2irU 2 



4tt 2 U 2 



(3.30) 



where 



O"ni<0 



-2i7r(t/(m+pi)- 



P2j 



< Pi < 1, C/ = Ui + i U 2 , U 2 > 0. (3.31) 



This is the result for the Casimir energy for T2/ZN, N > 2 with boundary conditions as 
in ( 3.20| ) and with p\^ taking the values required by the consistency conditions specific 
to each orbifold. Finally, V of ( 3.30| ) is an asymptotic series given by (see Appendix 
eq-O) 



V- 



e/2-2 



E 



2(-1)p 



■nl TT 



2p-2 



r 



P-2 + - 



C[-2p,p2-U 1 (n 1 +p 1 )} 



m 2 A 
(2ir) 2 U 2 



-K+pi) 2 



2-p- 



(3.32) 

where £[q, x] is the Hurwitz zeta function [26|, A = L 2 sin6, U2 = ImC7 = sin^, 9 = 2-k/N . 
This expression of T> is valid without any restrictions on the relative values of m, L or U. 

The quantity T> is of particular interest because it contains additional poles as e — > 0, 
and so potentially introduces new contributions to the UV divergent part of Vo in ( 3.3C| ). 



Note that if one of the sums (say that over ri2 > 0) in Vo of ( 3.29 ) were extended to the 
whole Z set of integers, the quantity T> given above would not arise due to the cancellation 
against the similar contribution to Vo, coming from 71-2 < 0. The latter would actually be 

P2 and Ui — > —U\. The sum of these 



1 



equal to V of (|3.32|) with the substitutions p2 
two contributions would then vanish 

i>(Pl,P2) +£>(Pl,l - P2, 



(3.33) 



4 with some exceptions in the case of T2/Z4, orbifolds, see later. 

5 For a general approach to computing traces on orbifold spaces see 
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since £[— 2m, x] + Q[—2m, 1 — x] = 0. In particular, this explains the absence of T> in 72 
and T2/Z2 where one of the KK sums was over the whole set Z. To conclude, the presence 
of T> in the potential Vq is due to the fact that both sums over the Kaluza-Klein modes in 



( 3.29 ) were restricted to positive/negative modes only. 



3.2.1 Ultraviolet Divergences 

Because the contribution D potentially introduces new UV divergences for orbifolds 
with N > 2, in this section we investigate their form in more detail in order to see what 
kinds of counterterms they require. In particular, we show that for T2/Z2N no new coun- 
terterms are required beyond those which already arise for the Z2 orbifold. To do so we 
consider in detail the case of I2/Z4, for which the analysis of T> is considerably simplified. 
For the remaining cases (with N = 3,6) the analysis follows the same technical 

steps as below, but is more involved and will be presented elsewhere |3(J . 

In the case of T2/Z4 which has U\ = 0, [7 2 = 1, the Hurwitz zeta function in ( 3.32j ) has 



no dependence on m, and this simplifies the identification of the additional poles. In the 
last bracket in eq.( |3.32| ) one can then use a binomial expansion (mi <C 1) or an asymptotic 
expansion (ml 3> 1) and following the technical details in Appendix |C|, eqs.( C-li| ), ( |C-12j ), 
(C-13) one obtains from ( 3.32j ) that, for T2/Z4 



~ A 2 fm 4 2m 2 (2vr) 2 1 (2vr) 4 ~ ^, ,1 

c 3 = C[-4,p2](l/2-/9 1 )+C[-4 )/9l ](l/2-p 2 )+2C[-2,p 1 ]C[-2,p 2 ], 

c 2 = (l/2- Pl )C[-2, P 2] + C[-2,Pi](l/2-p 2 ), d = (l/2- Pl )(l/2-p 2 ) (3.34) 

In the first expression T>f describes the terms 0(e ), and its exact value depends on whether 
mL is smaller or larger than unity, and is discussed later on. The zeta functions appearing 
in the coefficients are given by 

C[-2,x] = -\x{x-l)(2x-l) 

([-4, x ] = -— x(x - l)(2x - l)(3x 2 - 3x - 1) (3.35) 

oU 

Writing 

Vo = Vooc + Vg n (3.36) 
we therefore identify the following UV divergent terms: 

m 6 A m 4 a m 2 c 2 , vr 2 c 3 , . 

^ = 768^ + 16^ + lX7 + ^7 (3 - 37) 
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Notice that the structure of these divergences is valid independent of the relative size of 
m and A. For the Z4 orbifold we have seen that consistency of the boundary conditions 
requires we choose the value of p\ = P2 = p with p equal to or | , and so we must evaluate 
the coefficients with these choices. Since both ([— 2,p] and ([— 4, p] vanish when p = 
or p = \ , we see that for all such cases 

c 2 = c 3 = (3.38) 



leaving in Vo<x> only the divergences in the first and second terms in eq.( 3.37 ). 



The first term in ( |3.37 ) is a renormalization of the bulk cosmological constant and is 



present irrespective of the values of p\ or p2 ■ Its coefficient is 1 /4 the size of the similar result 
for the covering torus. Therefore, once this term is expressed in terms of the orbifold area, 
Ao = A/4, its coefficient is precisely the same as was found for T2 and T2/Z2, as expected. 
The term in fl3,37|) , proportional to c\, is a "brane" divergence, which renormalizes the 
brane tension. It is nonzero only for the case where p\ = P2 = 0, in which case c\ = 

To conclude, we find that the UV divergences for the Casimir energy due to compacti- 
fications on 72/^4 with discrete Wilson lines (pi,p2), have two kinds of divergences at one 
loop, similar to the case of T2/Z2. These have the form and coefficients required by the 
general heat-kernel analysis |22| and renormalize the bulk cosmological constant and the 
tension of branes localized at the orbifold fixed points. For the case of remaining orbifolds 
72/^3, ^2/^6, the analysis of the divergences of the quantity V of ( |3.32| ) is more involved 
since the Zeta function entering its definition will retain a n\ dependence. This makes the 
computation more tedious and the identification of the relevant counterterms more difficult 



to analyze in this case [30|. 



3.2.2 The finite part of Casimir Energy. 



For the finite part of the Casimir energy for the orbifold T^/^ one obtains in the limit of 



A <C 1 (see Appendix eqs.flCM^) and ([(Til)). 



2 A 2 

+ E 



2vr d 
45" 



__ p 2(l_ 5p 2_2pf + 6p?) 



(3.39) 



ni<0 



( , ^2 T • / n , 3 Fi +Pi\ 
(ni+pi) Li 3 (<7 n J + - 



2tt 



Ll4 (<7ni ) + ^2 Lis ) + CX - 



with the notation 

<J ni <o = e 27T ( ni+pl+ip2 \ 0<Pi<l (3.40) 
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Here Vf is an asymptotic series which for m 2 A <C 1 has the following expression (see 
Appendix 0, eq. flCM^) ) 

V f = vr 2 |2C[-2, pi]C[-2, p 2 ] HTire^- 1 ) + C[0, p 2 ] (c[-4, Pi] ln^re^i) + 2C'[-4, 1 - pi 



+ 



4C[-2, p 2 ]C'[-2, 1 - pi] + C[-4, p 2 ] (C[0, pi] ln(TrreT) + 2C'[0, 1 - p^ 



fc>3 



(-1)' 

fc! 



r[fc-2]C[-2fc,p 2 ]C[2fc-4,l -pi] L with r = (2tt) 2 /(p 2 ^)- (3.41) 



Since 6 for the orbifold T2/Z4 the Wilson lines have the values p\ =p 2 = 0, 1/2, Vf simplifies 
to give: 

3 



V 



p\=p-i 



=0 4"7T 2 



C[5], 



and 



V 



0. 



Pl=P2 = 



(3.42) 



Eqs.flpg), (|3l2D and also ( fDSf ), O), flOj ) give the final result for the Casimir energy 
for the orbifold T2/Z4 with discrete Wilson lines. Finally, with pi = p 2 = 0, 1/2 one has 
from eqs.(UD to fly|) 



P = 0,- 



(3.43) 



where y Ten is the result for the 2-torus given in eq.(2.10). Following closely these steps, 
one can also obtain from eqs.( |3.30| ), ( p. 32 ) similar results for T~2 /Zq and T2/Z3 orbifolds. 



3.2.3 Heavy- mass dependence. 

We discuss now the heavy mass dependence for the Casimir energy. For T2/Z4 with m> 1/L 



it turns out that the divergences in Vq are identical to those in eq.( 3.37 ). From Appendix 
0, eq- flOToD with @, ( |OT3l) , ( gig ) one obtains the full result for Vb for raL»l. Here 
we outline only the main behaviour which is 



V 



"1 4 Cl , . 2 2 1 m 6 .A 11 2/ 2 n 

ln(vre 7 2 m 2 /^) + l n (vr e m 2 /p 2 ) + • • • 



32tt* 



(3.44) 



with c\ = (1/2 — p) 2 and where the dots account for additional terms such as poly logarithms 
terms, identical to those in ( [3.39 ), and for (asymptotic series of) terms which are suppressed 
by inverse powers of m 2 A. The latter vanish in the special case of T2/Z4 with p = 0, 1/2, 
to leave only the (exponentially suppressed) polylogarithm contributions. 



4. Conclusions 

In this paper we compute the value of the Casimir energy for a very broad class of two- 
dimensional toroidal compactifications. These include the general case of T2 compactifica- 
tions with arbitrary boundary conditions for the 6D fields corresponding to the presence 



6 In (3.41) the derivative of Zeta function is taken wrt its first argument. 
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of arbitrary Wilson lines, as well as TiIZn orbifolds (also with Wilson lines) obtained by 
identifying points under Zjy rotations. Our calculations are explicit for a 6D scalar having 
an arbitrary 6D mass m, and we show how to extend these results to higher-spin fields for 
supersymmetric 6D theories. Particular attention was paid to regularization issues and to 
the identification of the divergences of the potential. The computation also investigated 
the dependence of the result on m, including limits for which m?A is larger or smaller than 
unity, (where A is the volume of the internal 2 dimensions). 

For the cases of T2 and T2/Z2, our calculation generalizes earlier results to include the 
dependence on an arbitrary complex structure, U, for the underlying torus. The potential 
obtained is likely to be useful for studies of the dynamics of these moduli, including their 
stabilization and their potential applications to cosmology ||, 21]. 

By carefully isolating the UV divergent part of V, we show that all of the divergences 
may be renormalized into a bulk cosmological constant (which gives a Casimir energy 
proportional to m e A) and - for the case of T2/Z2 - a cosmological constant (or brane 
tension) localized at the orbifold fixed points (which gives a Casimir energy proportional to 
m 4 ). Furthermore, these divergences agree with expectations based on general heat-kernel 
calculations, such as those recently performed for 6D compactifications in ref. 22]. For 
massive 6D scalar fields, m 2 A S> 1, the dependence on m of the finite part of the Casimir 
energy obtained in the modified minimal subtraction scheme, is exponentially suppressed. 

We present results for the Casimir energy also for T2/ZN orbifolds with N > 2, again 
including Wilson lines and any shape moduli which are allowed. The case T2/Z4 was stud- 
ied in particular detail. The UV divergences that emerge in this case again take the form 
required by general heat-kernel arguments, and can be absorbed into renormalizations of 
the bulk cosmological constant and brane tensions localized at the orbifold fixed points. 
The finite part of the Casimir energy was computed in detail and may be used for phe- 
nomenological applications. Finally, the technical tools of the Appendix can be used for 
other applications such as the one-loop corrections to the gauge couplings in gauge theories 
on orbifolds, in the presence of discrete Wilson lines. 
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Appendix 

A . Calculation of the vacuum energy in DR for 2D compactifications. 

We provide here details of the calculation of the vacuum energy. One has (d = 4 — e) 



v*( P1 , P2 )=^- d j2 



d d p 



In 



P 2 + Ml 



71 2 



-J£ \- 



(2tt)^ Jo ti+d/2 



/i is a finite, non-zero mass scale introduced by the DR scheme. A "prime" on a double 
sum excludes the (ni,n 2 ) = (0,0) mode. If a level (ni,n 2 ) is massless M niin2 = (for 
example if £ Z)), mathematical consistency requires one shift M% 



112 



by a finite non-zero m 2 = 5 p? (5 dimensionless) . This also helps us identify the scale (to) 
dependence of the divergences (poles in e). We use 



m: 



(2vr) 5 



,,^n 2 + p 2 -U(n 1 +p 1 )\ z ; U = U x + iU 2 = e w -4; .4 = L X L 2 sin (A-2) 

-AC/2 Li l 



The DR regularized sum in ( [A- 1| ) is re-written 



F*( Pl ,p 2 )=/i 4 C e £ / 



m. 2 ez ' 



-^*[A^ lin2 / M 2 +<5] 



(27T) 



4-e 



(A-3) 



with <5 = m / p . The calculation of is reduced to that of ^7* performed below. 



B . Series of Kaluza-Klein integrals and their DR regularization. 

We evaluate (in the text r -» (2ir) 2 /(p 2 A U 2 ), 5 -> m 2 / p 2 ) 

/ T£ e — *-l«2+P2-^(™i+Pi)l 2 e — 5 *, r>0 , U=U 1 +iU 2 , 5>0, pi GR(B-l) 

i7 £ * includes shape moduli effects (9^ir/2, L\^L 2 ), 5 shifts, and arbitrary "twists" pi wrt 
Ll,2- To evaluate J* one uses re-summation (|E-l| ); the integrand becomes 



2 ^ + i_„ rr„, |2 



e -irtT\n2+p2-U(n 1 +p 1 )\ 2 _ g -7r i r\n2+p2—Upi | _|_ g-?r t T|n2+P2-^(«.i+pi ) 

ni,n 2 GZ ri2GZ meZn2£Z 



E e -TTtT\n 2 +p2-U Pl \ 2 , 1 e -7rtrU"|(m- 

n 2 GZ v mez 

_i ^2 e -T^- 7r * rc/ i( n i+pi) 2 + 2 ™2(P2-{/i(pi+ni)) ^ B _ 2 ) 



mez n 2 ez 
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A prime on a double sum indicates that (ni, n 2 ) 7^ (0, 0) is excluded and a "prime" on a sin- 
gle sum excludes its n = mode. The three contributions above can be integrated termwise 
for any real pi, (given the presence of e - ^* 5 ). Accordingly, one has three contributions 

J*{p\,P2) = £i(pi,p 2 ) +/C 2 (/3i,P2) +/C 3 (pi,p 2 ) (B-3) 

defined/evaluated in the following: 
• Computing /Ci : 

1 poo jj 2 

Kx=Y,l -I^e-^I^-^IV^* = ?-(5 + T\p 2 -U Pl \ 2 ) 2 (B-4) 
+ ^(S + T\p 2 -U Pl \ 2 ) 2 ln 
+ t(5 + t U 2 p 2 ) Li 3 (e- 2 ^°)) + ¥1 (5 + r t/fp 2 )i Li 4 ( e - 2 ^ )) + ^ Li 5 ( e - 2 ^ ))+ C .c. 



(5 + r|p 2 - Upi 



|2\2 



8vr d 



(5 + rt/fp 2 



where "c.c." applies to the PolyLogarithm functions only. To evaluate K\ we first added 
and subtracted the n 2 = mode contribution. We then used a (Poisson) re-summation 
over n 2 , then the integral representation of modified Bessel functions ( |E-2| ) with (E-3), and 
the definition of the Polylogarithm hi a (x) flE-4| ) . Finally we used the notation 

1 



7 (0) = ^=(5 + rU 2 p 2 )2 - i(p 2 - U lPl ) 



(B-5) 



The divergence of IC± is that of the excluded n 2 = mode in JC\, which is in turn due to 
the absence of (ji\,ri2) = (0,0) in the definition of J* . 

• Computing JC 2 : We introduce the notation A pi =pi — [pi], 0<A pi <1, [pi]eZ. 
(a) . For < 5/ (rC/| ) < 1 we have 



dt 



-nt t C/f (ni+Pl) 2 —irSt 



1 



dt 



vWo * 7/2+e 

VrS/2+e 

— r[-V2- £ ] 



^ e -^t[«5 +r C/|(n 1+ A Pl ) 2 ] _ e -7rt(«5+rC7|p2) 
nieZ 



£ [ 5 + rC /| (ni + Api ) 2 ] 2 ^ _ (5 + rU 2 p 2 ) 2 

nieZ 



! + (<5 + t£/ 2 2 A 



pi' 



5 

7T 2 

^=r-5/2 



2„ 2\ 



+ 



(vrr[/ 2 )2- 
v 7 ^ 



£ 

fc>0 



r[fc-5/2-e] 



rC/f 



C [ 2 fc- 5 -2e, 1 +A P1 ] + (A Pl -A 



pi/ 



(B-6) 
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In the last step we used the binomial expansion 



T[k + s] 



£[a(n + C )2 + r s = a-*£^ 

n>0 fc>0 ' L J 



1 k 



C[2k + 2s,c], 0<q/a<l (B-7) 



Here a] with a / 0, —1, —2, • • • is the Hurwitz zeta function, (with £[q, a] = X^n>o( a + 
n)~ q for Re(g) > 1). Hurwitz zeta-function has one singularity (simple pole) at q = 1 and 
C[q, 1] = ([q] with £[</] the Riemann zeta function. The only divergence in /C 2 is due to its 
k = 3 term in (B-6), from the singularity of the Zeta function. In the remaining terms in 
the series one can safely set e = 0. Further 

1 



C[l-2e,l±A pi ] = 

r[l/2-e] = vr 1 / 2 (l + eln(4e 7 )) + 0(e 
x e = l + elnx + C(e) 



(B- 



we find for < 5/{tU%) < 1: 

7T 3 5 3 1 7T 3 5 3 



JC 2 



6r\U 2 \e + 6t\U 2 \ 
8vr 3 ,, , 



111 



4irTU2e 1+1p{Api)+ ^ { ~ Api) 



15^ 
-4vr 3 

r 5/2 



(5 + t£/ 2 A 2 ) 2 - (5 + t^ Pi ^)2 + (r^)5(C[-5, 1+ A„J + fl-5, 1- A P1 ]) 



5r 3 / 2 |C/ 2 | 3 (C[-3, 1 + A pi ] + C[-3, 1 - A p J) + vr 3 5 2 \U 2 \ (1/6 + A 



V? ^ (p + 3)! 



rf/f 



p+3 



C[2p+l,l + A Pl ] + C[2p + l,l- A-J 



(B-9) 



The divergence in /C 2 is due to n 2 = (i.e. Poisson re-summed zero mode wrt to the second 
dimension) in the presence of infinitely many KK modes of the first dimension (rii). It is 
thus an interplay effect of both compact dimensions. The condition of validity of the above 



result < 5/(tU£) < 1 gives 5 // < 1/Li 2 ■ If S <C 1 the result (|B-9|) simplifies considerably. 



(b). If <V(tC/ 2 )>1 or 5^>1 eq.( |B-7] ) does not converge. If so, JC 2 is reevaluated as below: 



r _ i v r * 



1 



E e 

nieZ 



-Trt t C/| (rai+pi) 2 —irSt 



-Trt^+rC/Km+pi) 2 ] _ -7rt(<5+rC/|p 2 ) 



1 



trC/ 2 c „ 

niSZ 



E 



^-+2i7rnipi — 7r<5t 



+ V j T T[-3 - e] (B-10) 
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where the last term originates in the hi = term of the series. We find (with (|E-2| 



/Co 



7T 3 <5 3 



1 



6t\U 2 \ Le 



+ In (vr^- 11 / 6 +4 x /?C/|<5l E 



2 (-3 ^ cos[27rnipi] 



nj>l 



8?r3 ,r rr2 2\5 



(B-ll) 

rapidly convergent if <5>rC/| or <53>1. This ends our calculation of /C 2 at large/small 5. 
• Computing /C3: 

; /"OO Jj 

"* " = -7rn|/(tr)-7rtr (7|(ni+pi) 2 +27rm 2 (p 2 -C/l(pl+ni))-7n5t ^g_-^2^ 



v nieZn 2 ez 



. g " "2/ 



Since /C3 is always exponentially suppressed at t —* and at i — > 00 it has no singularities. 
We can thus safely set e = 0. One finds 



z(n 1 )Li 3 (e- 2 ^" 1 )) + ^-(rz(n 1 )) 1 / 2 Li4(e- 2w ^ 1 ))+^Li 5 (e- 27r ^ 1 )) 



2(n x ) = (5 + rC/|(ni 



7 (m) = -^(S + rt/ 2 (ni + p^ 2 ) 1 / 2 - i{p 2 - U x {n x + p x )) 



+ c.c. 



(B-13) 



To evaluate /C3 we used the representation of Bessel functions eq.(|E-2|), then (E-5) and 
finally the poly logarithm definition in (|E-4| ). This result simplifies considerably if S -C 1. 



To conclude if 0<6/(tU$) < 1 we find for J* (with Q, Q), Q), Q, (^T|)) 

# = E 



ni,n2£Z 
_2 



e ~7rt T|ri2+p 2 -l/(rai+pi)| g" 71 "" 5 * 



£ (5 + r\p 2 - U Pl \ z Y + —{5 + r\p 2 - U Pl \ 2 ) 2 In 



vr 2 e 27 " 3 (J + t\ P2 - U Pl 



2\2 



+ T E 

niSZ 
87T 3 



z(m) Li 3 (e- 2 ^)) + JL( r *(m)) 1/2 U A {e~ 2 ^^) + ^ 14 (cr 2 **** 1 )) 
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! ^r 2 |f/ 2 | 5 [--A 2 1 (l-5A 2 1 -2A^ i ; 
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l-SOA^l + A^ 
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C[2p+l,l + A Pl ]+C[2p+l,l-A pi ] 



(B-14) 
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This restriction < 5/{tU 2 ) < 1 is required for the convergence of the calculation of K, 2 . 
The first line in J* is due to the absence of the mode (0,0). J* is well defined even for 
5 = if M nijn2 ^ 0. In such case the result is obtained by redoing the above calculation 
with 5 = or more easily, from the one above by formally setting 5 = 0. The above result 
simplifies considerably when 5 <C 1. 

Finally, if we have 8/{tU%) > 1, from eqs.flETTD, (gg). ( gg ), p^TT| ), ( p^l3| ) we find 



ni,n2SZ 
_2 



g -7ri T|n2+P2-^(ni+Pl)| 2 g- 71 " 5 * 



£ (5 + r|p 2 - IVilT + T (<5 + r|p 2 - l7pi|T In 



7T 2 e 27 ~ 3 (5 + r|p 2 -[/p : 



|2\2 



nieZ 



*(ni) LiaCe- 2 ^" 1 )) + A( TZ ( ni ))i/2 Li 4 (e- 2 ^ ni )) + % Li 5 (e~ 2 ^( ni )) 



+ C.C. 



^ r - +ln (vn^^j + V^f*t£ £E!^i^l * 3 ( 27 rn 1A /V(^f) } (B-lo, 



6r [/• 



ni>l 



and this concludes the evaluation of Using ( E-3 ), one shows that the last equation has 



the contributions from K% and from the polylogarithms suppressed if 5 3> rf/f an d ^ ^ T i 
to leave the first line and the term 7r 3 <5 3 / (6eT"| 1 ) as its leading behaviour for this region 
of the parameter space. 



• Adding to J* the effect of the mode (0, 0) the result is 

Je{pi,P2) = J*{Pl,P2) +Z e (pi,p 2 ) 

with 

Z € {fH,P2) 



(It 



o * 3+f 



e ~nt r\p 2 -Up\)\ 2 e -K5t 



-TT 



2e 



(5 + t\ P2 - U Pl \ 2 ) 2 - ^-(<5 + t\ P2 - U Pl \ 2 ) 2 In 



(B-16) 
(B-17) 

7r 2 e 2 7 -3 (5 + T | p2 _ c/w |2 ) 2 



J e is thus given by J* without the first line in QB-lQ or QB-15Q . 
Eqs. ( |2.8D , ( |2.10| ) in the text are then obtained from 

V(p 1 ,p 2 )=p 4 C e J_ e/2 (p 1 , P2 ) 1 with 5^m 2 /p 2 , t - (2vr) 2 /(^UU 2 ). (B-18) 
for the case i<C 1. 
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C . More series of Kaluza-Klein integrals for T 2 /Zn orbifolds. 
For general orbifolds T 2 /Zn one needs to evaluate (U = U% + i U2) 



E 



ni<0,n2>0 



dt 
1 



e -TrtT\n 2 +P2-U(n 1 +p 1 )\ 2 e -Kt5 



T[s}TT- s {r\n 2 + P2 -U{n 1 + Pl )\ 2 + 5) , r,S > 0, 0<pi )2 <l.(C-l) 

ni<0, rt2>0 



Eq.(p0|) in the text is then V a = -^/(2ir) A - e £_ i2 -±)(5^m 2 /n 2 , (2tt) 2 / (fi 2 AU 2 )) . 

To compute C s , the usual (Poisson) re-summation used in previous sections is not 
applicable given the restricted summation on n\,n 2 . The sum of the last line is actually 



a "truncated" Epstein function. To analyze C s , we follow the method in [32], for both 
non-zero pi 2 and complex U. For s = — (2 + e) this allows us to evaluate the Epstein 
function up to order 0(e 2 ), giving an expression for £_( 2 +e) U P to 0(e). We use that 



Ei[z; s; r, ci] = ^ [z + r(n 2 + ci] 

n 2 >0 



has the asymptotic expansion |l_5j] 



' m! 1 

m>0 



r[s + m] 



C[-2m,ci] + 



2i-s 



,1/2-s 



(C-2) 



T[g-l/2] 



27T S _£_! 
-I- T 2 4 7 2~ r 4 



^p s ~l cos(27reip) if^i (2^(z/r)5 J (C-3) 
p>l ^ ' 



Here £[q, a],a/0, —1, —2, • • • is the Hurwitz Zeta function, K s is the Bessel function (E-2). 
In the m = term in ( |C-3| ) one can use ("[0, ci] = 1/2 — ci even for ci = 0. One can use 
( C-3j ) recurrently for the 2D case |3^]. With the substitutions 

ci ->• ci(ni) = p 2 - Ui(n\ + pi); z -> z(m) = T[/f( ni + p^ 2 + 5 (C-4) 

in eq. flC-3| ) and after applying a summation over nj, one obtains from ( |C-3| ) £ s of (| 

£ s = B + C + P 



B = 7T~ 



2T[s- 1/2] 



E z ( n i)'" 



ni<0 



C = 2r a 4 ^ z(ni) 2 + 4p s 2 cos [27rci(ni)p] i^ s _i ( 27rp (z(rti)/r 
ni<0;p>l 

r[s + m] 



rti<0, m>0 



771 



L(-r) TO C[-2m,ci(ni)]z(ni) 



(C-5) 
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The series in T> is asymptotic [32]. To compute B one considers the cases S/(tU 2 ) < 1> 
> 1. In the following we take s = —2 — e. 

If < 5/(tJ7|) < I? one uses f° r ^ a binomial expansion of its term z(ni) 5 / 2+<E , as in 
eqs.( B-7p , ( |B-8| ) , and the comments thereafter to isolate the poles, to find 
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1 1 



+ 
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(p + 3)! 



rC/f 



p+3 



C[2p + l,l-pi] \ + 0(e) 



(C-6) 



If instead 5/{tU 2 ) — lj one uses f° r ^(ni) 5//2+e in £> the expansion eq. ( |C-3| ) , to find 
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(vr5) 5 / 2 ^ r[m - 5/2] /-tU : 
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r > ° 
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E 



C[-2m,l-pi] +0(e) (C-7) 



For C one uses the definition of Bessel functions K_ 5 / 2 an d of Li CT , to find 



C = ~J2 \<ni )Li 3 (e- 2 ^" 1 ) ) + 1- (rzim )) 3 Li^e" 2 ^™ 1 >) + |^Li 5 ( e ~ 2 ^ )) +c.c. 



ni<0 



4vr 2 ' 



+ ©(€) 



7 (m) = 4=0* + rC^Km + pi) 2 ) 1/2 - i(p2 - EMm + pi)) 



(C-8) 



Therefore if < <5/(tE7| ) < 1 one has from Q, 0, ( |08| ) 
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Eq.flOq) in the text is then V a = -^ 4 /(2vr) 4 - e £_ (2 _ f ){5^m 2 /p 2 , {2tt) 2 / {p 2 AU 2 )) 
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It remains to evaluate the series of V in dCl| ) in a form amenable to numerical evaluation. 
For this, its factor z(n\) under the sum can be expanded in ^(rCTf) if < <5/(rC/f) < 1 
by using the binomial expansion eq.( |B-7| ). If 8/(tU%) > 1 one uses instead the asymptotic 



expansion of eq. (|C-3|) . For our purposes s = — (2 + e), and then 



V = 7r 2+e V Y[m-2-e] y - f~C[-2m, P2~U 1 (n 1 + Pl )} 

^-^ ml 

ni<0, m>0 



5+TUl{ni+ Pl f 



2-m+e 



(C-10) 



Eq.(pg) in the text is V = 2A 2 ^/{2^~ e V\e^> -e/2, (2tt) 2 /(fi 2 U 2 A), 5^m 2 / P 2 j . 

• Case of 72/Z4 orbifold: In the following we restrict the calculation of T> to T2/Z4, when 
U% = 0, XJ<i = 1. If so, the argument of zeta function in T> does not have a n\ dependence 
and the sums over n\ and m can be easily performed. (For other orbifolds U\ 7^ further 
evaluation of T> is more tedious but very similar). 

(a). For < S/t < 1, after a binomial expansion flB£7| ) of last bracket in (|C^1C| ), 2? becomes 



771! 



2? = (vrr) 2+£ ^ r[m-2-e 
ni<0, m>0 

= (7rr) 2+e ^ ^j P ( 1)m +p r[p + m-2-6 



' C[-2m,p 2 ] J/T + K + pi) 2 



m>0,p>0 



ml pi 



■([-2m, pa] C[2p + 2m - 4 - 2e, 1 - pi] 



+ 



2e 



tt 2 5t 



C[-4, p 2 ] C[0, pi] + C[-4, pi] C[0, p 2 ] + 2C[-2, pi] C[-2, pa] 



} 2 5 2 
+ ^C[0, Pi] C[0, P2] + 2?/ + 0(e)(C-ll) 



used in flOg ), ( p7|) with P = 2/^ 4 L 4 /(2vr) 4 - e P(e^ -e/2, (2^) 2 /(^L) 2 , <5^m 2 //^ 2 ). 
2?/ in eq.( |r>dl] ) is the finite O(e ) part: 



vrr 



+ 



2C[-2, Pl ]C[-2, p 2 ] ln^-^+CfO, p 2 ] (C["4, pi] ln(™^-§) + 2C'[-4, 1 - Pl ]) 



4C[-2, p 2 ]C'[-2, 1 - pi] + C[-4, p 2 ] (C[0, pi] ln(vrTe^) +2C'[0, 1 - pj 



+ 7T 5t 



C[0, p 2 ] (c[-2, pi] ln^re^ 1 ) + 2C'[-2, 1 - p x ]) + C[-2, p 2 ] (([0, pi] ln^re ■ 



+ 2C'[0,1-Pi] 



p>0,m>0,p+m>3 



; , +m r[p + m 2] ^ + 2m _ 4 1 _ 2 

p! m! 
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This was used in ( gjgg) , ( PH) with = 2p 4 L 4 /(2^) 4 X>/(r -» (2^) 2 /(pL) 2 , <J -> m 2 /^ 2 ), 
after neglecting any mL <C 1 dependence. 



(b). In the case when S/t > 1 one uses in V of (C-IC) or the first line in ( C-ll| ), the 
asymptotic expansion eq. flC-3; ) . The results shows that the divergent part of V is identical 
to that in the last two lines in ( |C-11| ), while the value of Vf (O(e )) in ( |C-11| ) has now the 
expression 



V f = (nr) 2 Yl (-l) m+k Cl-2k,l- Pl ]([-2m,p 2 ] 



T\k + m-2\ /<J\2-m-fc 



m>0,fc>0,fc+m>3 



m\ k\ \t 



£_C[-2m >P2 ]LJ-^y r[m-5/2] 





m>0 


+ 


£ 




m>0 




7r 2 r 


+ 




2 



m! vr 



2r 2 (-vr) m /<hf-f 



m! vr 



C[-2m, p 2 ] ^ cospTrppi^-s ^ m _| (2vrp 
p>i 



C[-4, p 2 ]C[0, pi] + C[-4, pi]C[0, P2] + 2C[-2, pi]C[-2, p 2 ] 



ln(-7r5e 7 ) 



2r2 



11*6 



(C-13) 



+ ^ 2 r5(C[-2, Pl ]C[0, p 2 ] + C[-2, p 2 ] C[0, Pi]) ln^e^ 1 ) + —C[0, Pi] C[0, P2] ln(^" 2 ) 



To conclude, if <5/t > 1, the value of £_( 2+(E ) is given by 

£-(2+e) = S + C+P/ 



+ ^(C[-4, P2] C[0, Pi] + C[-4, Pi] C[0, P2] + 2C[-2, pi] C[-2, p 2 ] J 

+ ^ [C[0, Pi] C[-2, p 2 ] + C[-2, pi] C[0, pi\ I + ^C[0, Pi] C[0, P2] + O(e) (C-14) 

B is given in eq. (|C-7| ), C in eq. flog) , while Vf is that of eq.(^T|). 

Eq. (|C-14| ) concludes the calculation of £_( 2+e ) for T2/Z/1 for <5/r > 1. Eq. (|C-9| ) with 
dCTD , gives £_ {2+e) for < S/t < 1 again for T 2 /Z 4 . Eq.flCj) with @, Q, 

(jpj), give £_ (2+£) for any £/. 



With the expression (C-14) for £, eq.(3.44) in the text is then 

M 4 „ /. m 2 (2vr) 2 



r — > 



(C-15) 
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D . Orbifolds, Fixed points and discrete Wilson lines. 

The lattice of T 2 /Zn orbifolds is generated by (1,0 with £ = i for Z 2 , Z4 and £ = e 2l7r / 3 
for Z3 , ^6 . The group Zjv of discrete rotations has N elements Pjy , < n < N — 1 with 
Pjy = 1. Their fixed points are 

Z 2 : P 2 : z f . p . = 0, 1/2, £/2, (1 + OA £ = i 



Z3 : 


^3,^1: 




= 0, 


(2 + 0/3, (1 + 20/3, 




Z 4 : 




Z /.p. 


= 0, 


(1+0/2, 






p2 . 
M • 




= 0, 


(l + OA e/2, 1/2. 




Z 6 : 


-P6,P 6 5 : 




= 0. 








p2 p4 . 




= 0, 


(2 + OA (1 + 20A 






p3 . 




= 0, 


1/2, e/2, (1 + 0/2- 





_ p 2i7r/3 



£ = e 2ivr/3 



(D-l) 



The usual orbifold action (g) and that of Wilson lines (Ti^) are given by 

$»(z + l) = Ti$ 9 (z), 
$9( z + r ) = T 2 $ 9 (z) 

$ 9 (rz) = 5 $ 9 (z), with T = e 2i *' N 

One has that 

$f( z + T ) = ^9^ T -i z + !)) = 5Tl $9( T -i z ) = gT l g N ~ 1 <$> g {z) 
Using the definition of T 2 , then 

T 2 g = gTi 



(D-2) 



(D-3) 



(D-4) 



One can further assume that the orbifold action g and the Wilson lines Tj commute, then 
Ti = T 2 = T, and with T< = e 2i7rft one finds (modulo Z) that Pl = p 2 = p for any T 2 /Z N . 

The case of Z 3 orbifolds: (£ = e 2iw / 3 ) 

<5> 9 {t 2 (z + 1)) = 5 2 Ti$ 9 (z) 

$ 9 (r 2 (z + 1)) = Tf 1 $ 9 (r 2 (z + 1) + 1) = Tf 1 T 2 - 1 $ 9 (r 2 z) = Tf 1 ^" 1 g 2 <$> 9 {z) (D-5) 
A solution to this is 

g 2 T^^^-lgl (D _ 6) 
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or, assuming [g, T] =0 giving T\=Ti=T 

T 3 = l, => p = 0, 1/3, 2/3. (D-7) 

where T = exp(2m / o). Further, for the fixed points 

a), if z f/p . =0, =► $ 9 (0) = $ 9 (t0) = 5 $ 9 (0) 

6). i/ Z/ . p . = (2 + 0/3, Tg&{z f . p .) =&{z f . p .) 

c). if ^. = (1 + 20/3, T 2 5 $^z / . p .) = ^( Z/ . p .) = T 5 2 ^(z / . p .) (D-8) 

These are additional conditions which must be respected by the Wilson lines T, orbifold 
projections g and fields & 9 at the fixed points. The conditions can be respected by suitable 
relative choices for T, g, or trivially by requiring the fields vanish at these fixed points. 

The case of Z4 orbifolds: (£ = i) 

$ 9 {t 2 (z + 1)) = T^^ 9 (t 2 z + r 2 + 1) = Tf V $ 9 (z) 

$ s (r 2 (z + 1)) = 5 2 Ti $ s (z) (D-9) 

which gives 

g 2 T l =T{ 1 g 2 (D-10) 

or, assuming [3, Tj] =0 giving T\ = T2 = T one has 

T 2 = l, => p = 0,1/2. (D-ll) 

where T = exp(2i7rp). Further, for the fixed points 

a), if zf. p . =0, => $ 9 (0) = $ 9 (r0) = 5 $ 9 (0) 

6). i/ ^. = (1+0/2, Tg&(z f . p .) =&(z f . p .) 

c) . if z f . p . = 1/2, Tg 2 &(z f . p .) =T*f(-l/2) = <*> 9 (z/. P .) 

d) - */ z f . p .=S/2, => Tg 2 ^(z f . p .)=gT^(e/2) = ^(z f . p .)(D-l2) 

Similar to the Z3 case, the conditions can be respected by suitable choices for T, g, or 
trivially by requiring the fields vanish at these fixed points. 

The case of Z§ orbifolds: 

<5> 9 (t 3 (z + 1)) = g 3 T 1 ^ 9 {z) 

$ 9 (t 3 {z + 1)) = Tf 1 $ 9 (t 3 z + t 3 + 1) = Tf 1 £ 3 $ 9 (z) (D-13) 
since r 3 + 1 = 0. One solution is 

<7 3 T 1 = TfV (D-14) 
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Assuming [g, Tx] = which gives T\ = T 2 = T one has 

T 2 = 1, => p = 0,1/2. 

where T = exp(2i7rp). Further, 

§ 9 {t 2 {z + 1)) = g 2 T l ^{z) 

&{t 2 {z + 1)) = Tf 1 ^(r 2 z + r 2 + 1) = T^T 2 g 2 &{z) 
since r 2 — r + 1 = 0. One solution is 



g 2 T 1 =T 1 - 1 T 2 g 2 



With [5, Tj] = giving T\ = T 2 = T one has 



(D-15) 



T = 1 p = 0. 



(D-16) 



(D-17) 



(D-18) 



where T = exp(2i7rp). Thus, if [g,Ti] = 0, one concludes from ( p-15 ), ( D-18| ) that p = 0. 
Further relations at the fixed points exist, which can be found as in the case of ^2/^3,4. 

E . Mathematical Formulae and Conventions. 

We used the Poisson re-summation formula 



tA 1 h 2 +2inhcr 



A 

nez v ^ hez 

The integral representation of Bessel Function K u 



(E-l) 



dxx y - x e- hxP - ax 



-p 2 



2p 



Kv(2Vab), Re(b), Re{a) > (E-2) 



with 



K 3 (z » 1) = e 



1 3 3 

1+^T + ^T 



35 1 945 1 3465 1 
1 + ¥z + 128^2 + 1024^3 + 



The definition of PolyLogarithm 



Li a (x) = }2 — 



n>l 



(E-3) 



(E-4) 
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